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Gravitational radiation with roughly spherical wavefronts, produced by roughly spherical black 
holes or other astrophysical objects, is described by an approximation scheme. The first quasi- 
spherical approximation, describing radiation propagation on a background, is generalized to in- 
clude additional non-linear effects, due to the radiation itself. The gravitational radiation is locally 
T-H I defined and admits an energy tensor, satisfying all standard local energy conditions and entering 

. the truncated Einstein equations as an effective energy tensor. This second quasi-spherical approx- 

' imation thereby includes gravitational radiation reaction, such as the back-reaction on the black 

, hole. With respect to a canonical flow of time, the combined energy-momentum of the matter and 

• gravitational radiation is covariantly conserved. The corresponding Noether charge is a local gravi- 

D ' tational mass-energy. Energy conservation is formulated as a local first law relating the gradient of 

the gravitational mass to work and energy-supply terms, including the energy flux of the gravita- 
1^-^ ' tional radiation. Zeroth, first and second laws of black-hole dynamics are given, involving a dynamic 

surface gravity. Local gravitational-wave dynamics is described by a non-linear wave equation. In 
T-H . terms of a complex gravitational-radiation potential, the energy tensor has a scalar-field form and 

^ ' the wave equation is an Ernst equation, holding independently at each spherical angle. The strain 

CO to be measured by a distant detector is simply defined. 
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I , Gravitational waves and black holes are among the most popular and intensively investigated topics in physics at 
the turn of the millenium. Both are predictions of Einstein's General Relativity which involve essentially relativistic 
gravitational effects. Gravitational waves were originally introduced by Einstein himself Ql and the first black-hole 
solution was also almost immediately discovered though the term black hole was coined much later [Q. In the 
remainder of the twentieth century, astrophysical evidence eventually accumulated to the point where it is nowadays 
believed that black holes are not only common in the universe, but astrophysically dominant energy sources and 
' appreciable mass concentrations, being the final remnant of any sufficiently massive star, with supermassive black 
holes powering active galactic nuclei and lurking at the heart of most other galaxies, e.g. ||]. Gravitational waves are 
expected to make a similar transition from theory to observation with the operation of several new gravitational-wave 
detectors, e.g. Observation of gravitational waves from black holes would provide the first direct evidence for the 
existence of the latter, rather than the impressive but indirect evidence of their effect on other astrophysical bodies 
or surrounding matter. The combination of these two intertwined topics, gravitational waves and black holes, poses 
an exceptional challenge for theorists. 

For instance, an expected source of detectable gravitational waves is the inspiral and coalescence of a binary black- 
hole system. Although the earliest and latest stages are understood in terms of post-Newtonian ^ and close-limit 
approximations Q respectively, the coalescence is generally thought to be tractable only by full numerical simulations 
[g,|0|. Textbook theory simply does not suffice to understand the process in physical terms. Stationary black holes are 
well understood pl]-p^, but the black holes of interest are highly dynamical. Asymptotic or weak gravitational waves 
are well understood, respectively by Bondi-Penrose theory ]l^-pO| and Wheeler-Isaacson high-frequency (or merely 
linearized) theory pl|-p6[|, but strong gravitational waves produced by a distorted, rapidly evolving black hole are not 
understood at all. Indeed, it is sometimes argued that there is nothing to physically understand, merely complicated 
equations to numerically integrate, as gravitational waves cannot be localized The main purpose of this article is 
to address this lack of relevant theory by providing an astrophysically realistic approximation scheme in which both 
gravitational radiation and black holes are locally defined, along with their physical attributes, with each dynamically 
influencing the other. The approximation can be simply stated: it holds where the gravitational wavefronts and black 
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hole (or other astrophysical object) are roughly spherical. This does not contradict the classically quadrupolar nature 
of gravitational waves, as the wave amplitude at different angles may be arbitrarily aspherical. The wavefronts, not 
the waves themselves, should be roughly spherical. 

With such a general premise of rough sphericity, this quasi-spherical approximation scheme intuitively promises a 
wide range of validity, in particular including a just-coalesced black hole. It may also prove applicable to neutron 
stars or supernovas, or indeed any astrophysical situation which has rough spherical symmetry. Further, for any other 
process which can be enclosed by roughly spherical surfaces, it provides a mid-zone and, assuming isolation, far-zone 
approximation. Mathematically this was originally achieved ||27| by linearizing certain fields which would vanish in 
exact spherical symmetry, having made a decomposition of space-time adapted to the roughly spherical surfaces, 
henceforth called transverse surfaces. This will henceforth be called the first quasi- spherical approximation. Since 
there is no assumption of closeness to stationarity, the approximation holds for arbitrarily fast dynamical processes. 
The approximation has also been tested against the main expected source of asphericity, angular momentum, by 
applying it to Kerr black holes p8[ | : the error in the strain waveform is much lower than expected signals from binary 
black-hole coalescence. 

The approximation allows a local definition of gravitational radiation essentially because the decomposed fields 
naturally divide into quasi-spherical variables and wave variables, the latter satisfying a wave equation and yielding 
the Bondi news, or equivalently, the observable strain. Specifically, the gravitational radiation is encoded in the shear 
tensors of the outgoing and ingoing wavefronts. Then by retaining non- linear terms in these (previously linearized) 
fields, one would expect to obtain an approximation which is more accurate for the gravitational-radiation sector. 
This second quasi- spherical approximation is also presented here. Both first and second approximations share the 
remarkable feature that, to compute the observable waveforms, no transverse derivatives need be considered. The 
truncated equations form an effectively two-dimensional system, already written in characteristic form by virtue of the 
first-order dual-null formulation, to be integrated independently at each angle of the sphere. Numerically this is much 
easier to implement and computationally inexpensive when compared to the full Einstein system. Numerical codes 
exist for both first and second approximations |^^. Moreover, the wave equation can be written as an Ernst equation in 
terms of a gravitational-radiation potential. Thus the gravitational-wave dynamics is amenable to analytical methods. 

The first approximation describes gravitational-wave propagation on a background, though the background is neither 
necessarily spherical nor fixed in advance; this merely means that the quasi-spherical equations decouple from the 
gravitational-wave equation. There is no such decoupling in the second approximation, and therefore no background 
which is independent of the waves. Back-reaction of the waves on the geometry is thereby included. Specifically, one 
may define a gravitational-radiation energy tensor | p9[ | which acts just like a matter energy tensor in the truncated 
Einstein equations. Thus there is a fully relativistic inclusion of gravitational radiation reaction for dynamic black 
holes. For instance, if a black hole emits gravitational radiation, it backscatters to produce ingoing radiation which is 
absorbed by the black hole, thereby increasing its mass and area. This last property follows from a local second law 
@|l|,part of a general theory of black-hole dynamics p2) , where a black hole is defined by a type of trapping horizon. 
This local theory of dynamic black holes, not to be confused with textbook black-hole statics and asymptotics, is 
extended here by deriving quasi-spherical generalizations of the spherically symmetric first p3| and zeroth ||3^ laws 
of black-hole dynamics, involving local definitions of mass and surface gravity. 

The article is organized as follows. §11 reviews the formalism of dual-null dynamics, which provides a geometrical 
description of the wavefronts and the corresponding decomposition of Einstein gravity. §111 describes both first and 
second quasi-spherical approximations and the resulting truncated field equations. §IV locally defines the gravitational 
radiation and its energy tensor Q. A truncated Einstein tensor illustrates the role of as an effective matter energy 
tensor in the second approximation. §V locally defines dynamic black holes, along with mass m and surface gravity k. 
§VI derives an energy conservation law for the combined energy- momentum of the matter and gravitational radiation. 
This is also formulated as a first law relating the gradient of m to work and energy-supply terms, including the 
energy flux of the gravitational radiation. §VII gives zeroth, first and second laws of black-hole dynamics and various 
inequalities. §VIII introduces the complex gravitational-radiation potential, in terms of which Q takes a scalar-field 
form and the gravitational-wave equation takes an Ernst form. §IX defines conformally rescaled fields more suited to 
asymptotics, including a localized Bondi flux and a conformal strain tensor. §X concludes. 

II. DUAL-NULL DYNAMICS 

The wavefronts of outgoing and ingoing gravitational radiation form two families of null hypersurfaces, intersecting 
in the two-parameter family of transverse spatial surfaces. This geometry is described by the formalism of dual-null 
dynamics p^ , ^ , summarized in this section. Denoting the space-time metric by g and labelling the null hypersurfaces 
by , the normal 1-forms = —dx^ therefore satisfy 
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g-i(„±,n±) = 0. (1) 

The relative normalization of the null normals may be encoded in a function / defined by 

= -g^\n+,n-). (2) 

Then the induced metric on the transverse surfaces, the spatial surfaces of intersection, is found to be 

h = g + 2e^^n^ ® n~ (3) 

where ® denotes the symmetric tensor product. The dynamics is generated by two commuting evolution vectors u±: 

K,"-]=0 (4) 

where the brackets denote the Lie bracket or commutator. Thus there is an integrable evolution space spanned by 
The evolution derivatives, to be discretized in a numerical code, are the projected Lie derivatives 

A± - (5) 

where _L indicates projection by h and L denotes the Lie derivative. There are two shift vectors 

s±=_Lm±. (6) 

In a coordinate basis u_, ei, 62) such that u± — d/dx^ , where = d/dx°' is a basis for the transverse surfaces, 
the metric takes the form 

g = hab{dx"- + s'ldx+ + s" da;") (g) {dx'' + s+dx+ + stdx^) - 2e^^dx+ (g) dx~ . (7) 

Then (/i, /, s±) are configuration fields and the independent momentum fields are found to be linear combinations of 

6»± = 5i±*l (8) 

a±^ ±L±h-e±h (9) 

= L±f (10) 

hi[l-,l+]) (11) 



where * is the Hodge operator of h and L± is shorthand for the Lie derivative along the null normal vectors 

l±^u±-s±^e-fg~^{n'f) (12) 

which will be assumed to be future-null. Then the functions 9± are the expansions, the traceless bilinear forms a± are 
the shears, the 1-form to is the twist, measuring the lack of integrability of the normal space, and the functions are 
the inaffinities, measuring the failure of the null normals to be affine. The fields {d±, (j±, i'±,uj) encode the extrinsic 
curvature of the dual-null foliation. These extrinsic fields are unique up to duality ± 1— s- =p and diffeomorphisms which 
relabel the null hypersurfaces, i.e. dx^ i— > e'^^dx^ for functions A±(a;^). 

The dual-null Hamilton equations and integrability conditions for vacuum Einstein gravity have been given previ- 
ously in a different notation . They are linear combinations of the vacuum Einstein equation and a first integral 
of the contracted Bianchi identity. It is straightforward to add matter: denoting the projections of the energy tensor 
T by T±± = T{l±,l±), r+_ = T{1+, /_) and T± = ±r(/±), the resulting field equations are 

L±0± = -i^±e±-^9l~l\a±\^~87rT±± (13) 

L±0^ = -6+9^ - e-f (ii? - \^Df ±lu\^ + D- {^Df ± w)») + 8nT+^ (14) 

L±iy^ = i(T+ : al - \9+9^ - e'f {\R - 3|w|2 + ijD/p ± t^tJ . of) + 8n (T+_ + ie-^/i" : T) (15) 

±L±LJ = -9±Lo ± i(D« • CT± - D9± + Dv± - 9±Df) =F 87rT± (16) 
_LL±cr^ =a+-h^ ■<7^± \{9+cT^ - 9^cT+) + 2e~f {{\Df ±uj)® {\Df ±uj)-D® {^Df ± uj)) 

~e-f {\\Df ±uj\'^-D- {\Df ± LuY) h + 8TTe-f (_LT - \{h'^ : T)h) (17) 

_L (L+.s_ - i-s+) = 2e-fh^ ■ uj + [s-,.s+] (18) 

±L±h = 9±h + a± (19) 

L±f = y± (20) 
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where D is the covariant derivative and R the Ricci scalar of h, a dot denotes symmetric contraction, a colon denotes 
double symmetric contraction, a sharp (ft) denotes the contravariant dual with respect to = (index raising), 
a flat (b) will denote the covariant dual with respect to h (index lowering), |wp = lu ■ u)^ and \a\^ — a : crK Units 
are such that Newton's gravitational constant is unity. This is the Einstein system in first-order dual-null form; the 
equations will simply be called the field equations. The second-order version obtained by eliminating the momentum 
fields and symmetrizing in are the Einstein equations themselves. 



III. QUASI-SPHERICAL APPROXIMATIONS 



The dual-null fields and operators fall into two classes, those which vanish in spherical symmetry, {(J±,lu, s±, D), 
and those which generally do not, {d±, i'±, h, f, A±) [ ^J36| . The first quasi- spherical approximation therefore 
consisted of linearizing in (ct±, w, s±, D). It was then noticed that the gravitational radiation is encoded in the shear 
tensors a± |29), as will be explained in more detail in the following. This suggests that retaining non- linear terms in 
a± would give a more accurate approximation for the gravitational-radiation sector of the theory. Thus the second 
quasi- spherical approximation consists of linearizing in (w, s±, D) only. For a given matter model, one would also have 
to decide which decomposed matter fields to linearize. In this article, no specific matter model will be fixed, but the 
matter energy tensor T will be retained for generality and assumed to be consistently truncated. 

More formally, one may introduce two expansion parameters into the full field equations, ei preceding a± and eg 
preceding {u!,s±,D). Then both approximations ignore terms o(eo), whereas the first approximation ignores terms 
o(ei). Then ei measures the strength of the gravitational radiation, whereas eo measures other asphericities, due to 
angular momentum or other transverse effects. The terminology "first and second order" has been carefully avoided 
because the true second-order quasi-spherical approximation would be full Einstein gravity, since the only non-linear 
terms in the full field equations are second-order in the dynamical fields and operators. However, one may say that 
the second approximation is second-order in the gravitational radiation. 

It is useful to decompose the transverse metric h into a conformal factor r and a transverse conformal metric k by 

h^r'^k (21) 

such that 



A±*l = (22) 

where * is the Hodge operator of fc, satisfying *1 ir"^ . The Ricci scalar of h is found to be 

R = 2r-^{1- D'^lvLr) (23) 

by using the coordinate freedom on a given surface to fix k as the metric of a unit sphere. One may take quasi-spherical 
coordinates x"' — ip) on the transverse surfaces to obtain the standard area form of a unit sphere: 

il = sm-ffd-d Adip (24) 

where A denotes the exterior product of forms. Then r is the quasi- spherical radius. Approximations for rough 
cylindrical or plane symmetry could similarly be produced. A useful truncation identity, holding in both first and 
second approximations, is 

A± = -LL±. (25) 

This will be used throughout the article without further reference, with A± rather than L± appearing explicitly. 

In the first approximation, the truncated field equations decouple into a three-level hierarchy. The first-level 
equations 

(26) 
(27) 

^el - 8^r±± (28) 

+ 87rT+_ (29) 
B7r(T+_ + ^e-f'h^ : T) (30) 







A±/- 


>^± 


A±e± = 




A±0^ = 




A±v^ - 


-^0+9. 
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have the same form as in spherical symmetry. These equations, the quasi-spherical equations, therefore determine a 
quasi-spherical background. The second-level equations 

A±fc = r"V± (31) 
AiCTzp = 5(6'±cr^ - 0^a±) + cr+ • • cr_ + 87re-/(_Lr - : T)h) (32) 

constitute a wave equation for fc, describing the gravitational-wave propagation, as explained in detail in the following. 
The quadratic shear term in the shear propagation equation (^2|) was omitted in the original reference [ p7| . The third- 
level equations for (w, s±) need not be solved for the radiation problem. This is because, fixing to be the outgoing 
direction, the Bondi news at null infinity 3^ is essentially (Tzp/r [ ^Tt - This determines the strain to be measured by 
a gravitational-wave detector, as explained in the penultimate section. 

In the second approximation, the truncated field equations also decouple, this time into only two levels, with the 
last level for {uj, s±) again being irrelevant to the radiation problem. The remaining equations are 



A±r = ir6'± (33) 

A±/ = (34) 

A^fc = r-^cr± (35) 

A±0± - -iy±e± - \el ~ i|a±|2 - 87rT±± (36) 

A±6'^ = - e-^'r-^ + 87rT+_ (37) 

A±i/^ = -^0+9- - e-fr^^ + \a+ : at + 87r(T+_ + \e-fh^ : T) (38) 

A±cr=F = \{0±c^^ - Ot^±) + cr+ • ■ + 87re-/(_Lr - i(/i« : T)h). (39) 



The dual-null initial-data formulation is based on a spatial surface S and the null hypersurfaces S± locally generated 
from S in the u± directions. The structure of the truncated field equations shows that one may specify (9±, r, /, k) on 
S, ((T-|_,i^+) on S+ and (cr_,i^„) on S_. In particular, the initial data is freely specifiable. In summary, the vacuum 
system consists of nine first-order differential equations and their duals. 

Mathematically, the difference between first and second approximations is that in the first approximation, the 
equations for the quasi-spherical variables {9± ,i'±,r, f) decouple from the equations for the wave variables {a± , k) . 
Physically this describes gravitational-wave propagation on a quasi-spherical background. The background is not 
fixed in advance and need not be spherically symmetric, so even the first approximation is widely applicable. There is 
no such decoupling in the second approximation: the gravitational- wave terms {a± , k) now enter the equations for the 
quasi-spherical part of the geometry, so there is no longer a background which is independent of the waves. Physically 
this corresponds to including radiation reaction, as clarified in the next section. 

Nevertheless, both first and second approximations share the remarkable feature that, to compute the observ- 
able waveforms, no transverse D derivatives need be considered. The truncated equations form an effectively two- 
dimensional system, to be integrated independently at each angle of the sphere. Physically this means that the 
observed gravitational- wave signal depends only on the line of sight to the source, surely a plausible result. Moreover, 
by virtue of the dual-null formulation, the equations are already written in characteristic form, the mathematically 
standard form for analysis of hyperbolic equations. Numerical implementation is consequently straightforward and 
computationally inexpensive. Numerical codes exist for both first and second approximations pq|. 



IV. GRAVITATIONAL RADIATION: LOCAL ENERGY TENSOR 



The truncated equations in the second quasi-spherical approximation, (33)-(|39|), differ from those of the first, ( |2(: 
(p2|), by terms quadratic in the shear tensors a±, which appear additively with terms involving the energy tensor of 
the matter. Specifically, the second approximation may be obtained from the first by replacing the matter energy 
tensor T with T + O, where 

e±± = (7± : CT^/327r (40) 

e+_ = (41) 

_Le = ef {a+ : at)h/ 3277. (42) 

Thus 8 playes the role of an effective energy tensor for the gravitational radiation. Written covariantly in terms of 
the transverse conformal metric k, this defines the energy tensor of the gravitational radiation: 
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fc'a/a = (,4d) 

oZtt 

where {a, f3) = k"'''k'''^aacPbd is the transverse conformal inner product and A = dx^A^ + dx^ A^ is the l-form 
evolution derivative. 

More formaUy, in terms of the expansion parameters eo and ei, one may introduce a truncated Einstein tensor C 
defined in terms of the fuU Einstein tensor G by 

C±± = hm lim G±± = -A±0± - v±e± - ^91 (44) 

C+_ = hm hm G'+_ = A+e_ + 9+9 _ + e^-^r"^ (45) 

ei— >0 eo— »0 

: C = hm hm : G = {2A+v^ - 2A+9^ - 9+9^) (46) 

£1^0 eo^O 

(7= hm G' = e-^(A(+(7_)-(T+-/i»-CT^) (47) 

where C = J-G — ^{h^ ■ G)h and the G± components are irrelevant to the radiation problem. Then 

G = SttT (48) 
are the truncated Einstein equations of the first approximation, obtained from ([2^)-(|3^), whereas 

G = 87r(r + e) (49) 



are the truncated Einstein equations of the second approximation, obtained from (p3[)-(|39|). This demonstrates that 
is an energy tensor for the gravitational radiation, in the sense that it is included as an effective matter energy tensor 
in the second approximation. Gravitational radiation reaction, the back-reaction of the radiation on the space-time, 
has thereby been included. 

It should be stressed that (i) mathematically, O is a genuine tensor, but depends on the dual- null foliation, not 
just on the space-time; (ii) the physical interpretation of O as energy requires the quasi-spherical approximation to 
be valid, meaning that the transverse surfaces must indeed be roughly spherical. This will not be made precise here, 
as the range of validity of the approximation is not clear in advance and best explored in applications. The intuitive 
meaning of roughly spherical should be clear by any standards. 

In summary, the quasi-spherical approximation allows a local definition of the energy-momentum-stress of gravita- 
tional radiation, and therefore of the radiation itself: gravitational radiation is present at a given point if and only if 
G is non-zero there. With the orientation such that u+ is the outgoing null direction, there is outgoing radiation if 
and only if A^k (equivalently (t_) is non-zero, and ingoing radiation if and only if A+fc (equivalently (t+) is non-zero. 
The terminology gravitational radiation rather than wave is generally preferable, since A±k need not be oscillatory. 
Instead, frequency spectra for ingoing and outgoing radiation may be defined by Fourier transformations to frequency 
/+ and /_ respectively: 

fc^(/±) = / e^2^*^±^*fc(a;±)dx± (50) 

where 7 is a curve of constant (a;^,i?, <p). If the Fourier transform is peaked in frequency space, one may say that 
there is a gravitational wave. In contrast, gravitational radiation is defined even when there is no typical frequency. 
This clearly indicates that the approximation has a different physical basis to that of the Isaacson high-frequency 
approximation p^ , p5[ |, which is usually quoted to make sense of linearized gravitational waves. 
The non-zero components of Q may be written as 

e±± = ||A±fc||V327r (51) 
_Le = e-''(A+fc,A_fc)/i/327r (52) 

where | \a\ p = (a, a) is the transverse conformal norm. Then the Q±± components have a similar form to those of the 
gravitational- wave energy tensor of the high-frequency approximation p^ , p5[ | , with k replacing the transverse traceless 
metric perturbation. However, the high-frequency approximation requires averaging over several wavelengths and has 
no term proportional to g. In this connection, it seems that the quasi-spherical situation allows a natural choice of 
transverse traceless gauge and eliminates the need for averaging. Earlier attempts to construct pseudotensors for 
gravitational waves by Einstein and others might be converted to genuine tensors by similar gauge-fixing adapted to 
the transverse surfaces, but currently such pseudotensors are generally not accepted. 
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Apart from the fact that a gravitational wave generally has two polarizations, as encoded in the two independent 
components of fc, 8 is analogous to the energy tensor of Einstein- Rosen gravitational waves jss), which takes the 
massless Klein-Gordon form in terms of a gravitational potential generalizing the Newtonian potential. In both 
cases, there is generally a transverse radiation pressure, given by -L0, produced by a combination of ingoing and 

outgoing gravitational radiation, as well as the expected radial radiation pressure. The fact that = may be 

interpreted as meaning that gravitational radiation is purely radiative and workless, as for the massless (but not 
massive) Klein-Gordon field. 

As an energy tensor, Q satisfies the strong, dominant, weak and null energy conditions [^9| , as follows. The dominant 
energy condition, requiring that an observer measures future-causal momentum, may be stated as T(a, (3) > for two 
future-causal vectors a, 13. Writing the vectors as 



a 



^a+l++a^l^ + {2e~fa^a-Y/'^a _La = a (53) 
(3 = b+l+ + b-l^ + {2e-^'b+b-f'% l.b = b (54) 

the future-causal conditions imply > 0, b^ > 0, \a\ < 1, |6| < 1 and therefore < a ■ b'^ < 1. Then 
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327re(a,/3) =a-6^ {a+b+)^/^A+k + ia-b-)^^^A^k + {1 - a ■ b^) {a+b+\\A+k\f + a-b-\\A^k\f) > 0. (55) 

Thus Q satisfies the dominant energy condition and therefore also the weak and null energy conditions. The strong 
energy condition states that T(a, a) > where T = T — ^{g~^ ■ T)g. Since 

e„ . (^4^) ,56, 

this follows more simply as 

327re(a, /3) = | |a+A+fc + a" A_fc| |^ > 0. (57) 
Thus by all standard measures, gravitational radiation carries positive energy^ a physically expected property. 



V. BLACK HOLES: MASS, SURFACE GRAVITY 

A general definition of dynamic black holes was proposed previously [^,0. A transverse surface is said to be 
trapped if 9^9^ > 0, marginal if 0+ = or 0_ = and untrapped if 9^9^ < 0, also previously called mean convex [ p7[ . 
Such definitions may also be applied independently at each spherical angle. The surface is future trapped if 9± < 
and past trapped if 9± > 0, or future marginal if 0_ < (for = 0) and past marginal if 9- > 0. An untrapped 
surface has a preferred spatial orientation ||37| : ii 9+ > and 9- < 0, then A+ is outward and A_ inward. Similarly 
any spatial direction normal to the surface is outward or inward as its component along A+ is positive or negative 
respectively. Future or past marginal surfaces similarly have a preferred spatial orientation. A trapping horizon is a 
hypersurface foliated by marginal surfaces. It is an outer trapping horizon if < (for 9^ — 0) and an inner 

trapping horizon if L_9^ > 0, where the dual- null foliation is adapted to the marginal surfaces. A future (respectively 
past) outer trapping horizon is proposed as the local, dynamical definition of a non-degenerate black (respectively 
white) hole. A degenerate black hole is one for which 9+ is decreasing in the L_ direction, but L^9+ is not strictly 
negative. 

These ideas of gravitational trapping may also be expressed in terms of the quasi-spherical radius r, as in spherical 
symmetry |3^ , as follows. A transverse surface is trapped, marginal or untrapped as A^r is temporal, null or spatial 
respectively. If AV is future (respectively past) temporal or null, then the surface is future (respectively past) trapped 
or marginal. On an untrapped or marginal surface, an achronal (spatial or null) normal direction is outward or inward 
as r is increasing or decreasing respectively. A trapping horizon is outer or inner as A^r > or A^r < respectively, 
where A^ = A" • A. 

In spherical symmetry, there is an active gravitational mass-energy with many desired properties |^,^ . A quasi- 
spherical mass-energy may be defined similarly as 

m = ir(l - Ar • AV). (58) 

This is the simplest generalization which has the property that a transverse surface is trapped, marginal or untrapped 
if and only if r < 2m, r — 2m or r > 2m respectively. Physically, m includes the energy of the gravitational field, as 
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will be seen below in the first law. Similarly, the spherically symmetric m yields, in a post-Newtonian expansion, the 
Newtonian mass, gravitational potential energy, kinetic energy and thermal energy of a thermodynamic material [^ . 
A canonical flow of time is generated by the vector 

^ = (*dr)« = (*Ar)» (59) 

where * is the Hodge operator of the evolution space, *1 = dx^ A dx^ . One may define ^ equivalently, up to sign, 
by 

C • Ar = (60) 
C • = -Ar • AV (61) 
_LC^ = 0. (62) 

Then a sphere is trapped, marginal or untrapped as ^ is spatial, null or temporal respectively. In particular, trapping 
horizons may equivalently be defined as hypersurfaces where ^ is null. In these and other ways, ^ is analogous to the 
Killing vector of a stationary space-time or the Kodama vector of a spherically symmetric space-time ]3^ , |3^ . 
Denoting the Hodge operator of g by *, there is the truncation identity 

* 1 = *1 A *1 = *1 A ir^ (63) 

since the terms in the shifts s± are quadratic. Then 

V • C = = i^d{dr A Sr^) = r-^^d{r^dr) = 0. (64) 

Thus ^ is covariantly conserved. This Nocther current therefore admits a Noether charge 2-form such that 

^^^*dQ^. (65) 

Reversing the above calculation yields 

= 4r3. (66) 
Thus the integrated charge over a transverse surface is 

y = j>Q^^\§*r^ (67) 

which is the average radial volume. 



A dynamic surface gravity k may be defined, generalizing that of the spherically symmetric case 33 1, by 

T(e • d^^) = KAr (68) 

where T indicates projection onto the evolution space. That is, the 1-forms on each side of the equation are pro- 
portional in the quasi-spherical truncation, and k is defined as the proportionality constant. Since — ±Ar on a 
trapping horizon, this reduces on such a horizon to 

T(C • dt) = ±Kt (69) 

which is analogous to the usual definition of stationary surface gravity, with ^ replacing the stationary Killing vector. 
A calculation shows that 

K, = iA^r. (70) 

Therefore a trapping horizon is outer or inner as k > or k < respectively. Also, k vanishes where a black hole 
becomes degenerate, a desired property of surface gravity. Note that k is defined everywhere in the space-time, not 
just on horizons. Physically, k is a gravitational acceleration, a relativistic version of the Newtonian gravitational 
acceleration or force per unit mass. 
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VI. ENERGY CONSERVATION: FIRST LAW 



The truncated field equations (|33|)-(p9|) imply that the mass-energy m is propagated as 

A±m = inr^ef (T+„A±r - {T±± + e±±)A^r) (71) 

where the terms in Q ( ^ ) appear in the second but not the first approximation. This may be written more covariantly 
as the unified first law 

Am = 47rr^(V' + wAr) (72) 
in terms of the work density (an energy density) 

w = e^T+_ (73) 

and the energy flux (a 1-form) 

V; = V4T]+V[e] (74) 

which has been divided into contributions from the matter and the gravitational radiation: 

ij[T]=T{T- AK)+wAr (75) 
ip[e] = T(e • A«r). (76) 

The terms in the first law (fz^ ) involving -0 and w may be interpreted as energy-supply and work terms respectively, 
analogous to the heat supply and work in the classical first law of thermodynamics pO| . The unified first law was 
so called because it was shown in spherical symmetry that projecting it along the flow of a thermodynamic fluid 
yields a first law of relativistic thermodynamics, while projecting it along a trapping horizon yields a first law of 
black-hole dynamics It also includes the Bondi energy-loss equation at ^J"*", with m reducing to the Bondi mass, 
r^V reducing to the Bondi fiux and w reducing to zero. The quasi-spherical first law ( |7^ ) is unified in another sense: 
it includes the energy flux of the gravitational radiation, as also occurs for the first law in cylindrical symmetry . 
The energy-momentum density of the matter, referred to the canonical flow of time, is the vector j[T] given by 

/[T] = -T(r.O. (77) 

In spherical symmetry, the analogous j [T] is also conserved |^,^ . This also holds in the first but not second quasi- 
spherical approximation, the physical reason being that gravitational radiation carries energy. The energy-momentum 
density of the gravitational radiation is the vector given by 

/[e] = -T(e.e). (78) 

The energy fiuxes are essentially duals of the energy-momentum densities: 

*j^[T] = ?/;[r] + wAr (79) 

*j^[e] = 0[e]. (80) 

It follows from the first law that the combined energy-momentum density 

j=j[T]+m (81) 

is given by 

. ^ {^dmf ^ (*Am)« ^^2^ 

Then, again using the truncation identity (|63|), j is covariantly conserved: 

„ , u -kdidm A il) *ddm , , 

Physically this represents conservation of energy for the gravitational radiation and matter combined. The Noether 
charge 2-form Qj associated with the Noether current j is given by 
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j'=*dQ,. (84) 

Reversing the above calculation yields 

- ^ (85) 
or essentially the mass-energy m. The integrated charge over a transverse surface is 

m = <j) Qj — ^ j) *m (86) 

which is the average mass-energy. Thus the first law ( [72| ) is a first integral of the energy conservation equation ( p3[ ) . A 
general method for defining mass-energy by Noether charges has been recently suggested pl| . In summary, there is an 
energy conservation law for the combined energy-momentum of the matter and gravitational radiation. Physically this 
shows how energy may be locally transferred between matter and gravitational radiation, the balance being accounted 
by the active gravitational mass-energy m. 



VII. BLACK-HOLE DYNAMICS: LAWS, INEQUALITIES 

A general local theory of black- hole dynamics was initiated a few years ago and recently reviewed . The 
second law of black-hole dynamics ]3C| ] states that the area element of a future outer trapping horizon is non-decreasing, 
assuming the null energy condition. Thus 

r' > (87) 

where the prime denotes the derivative along a vector generating the marginal surfaces, with orientation such that, 
when r' — 0, it is future- null. This is related to the signature law that an outer trapping horizon is achronal 
(spatial or null), being null if and only if r' = 0. The local second law implies an integral version: the area A = f ir^ 
of the marginal surfaces is also non-decreasing, A' > 0. Since r = 2m on a trapping horizon, the second law also 
implies 

m' > (88) 

so that the black-hole mass is also non-decreasing. A similar monotonicity property of m holds in an untrapped 
region, as follows. It is straightforward to show that ip[Q] is past (respectively future) causal in future (respectively 
past) trapped regions, and outward achronal in untrapped regions. The same holds for assuming the null energy 
condition. Also, the dominant energy condition implies w > 0. Therefore projecting the first law ( [72|) along an 
outward achronal direction, m is non-decreasing in that direction. This generalizes the monotonicity property of m 
in spherical symmetry p6| and correspondingly yields two important inequalities. Firstly, positivity: in an untrapped 
region achronally outward from a regular centre, 

m > (89) 

since m vanishes at a regular centre. Secondly, in an untrapped region achronally outward from a marginal surface 
with r = rg, 

m > ro/2 (90) 

which is a local version of the Penrose inequality for black holes. These inequalities extend to the asymptotic mass 
at 3^ or spatial infinity i^ in an asymptotically flat space-time, as in spherical symmetry js^. 

An explicit expression for the surface gravity ( [70[ ) follows from the truncated equations of either first or second 
approximations, (|3^-(^): 

m 

K = — — ■iirrw. (91) 



In vacuo, this has the same form as Newton's inverse-square law of gravitation, combined with local equivalence of 
inertial and passive gravitational mass. The dominant energy condition yields the inequality 
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m > r'^K. (92) 
Also, in an untrapped region or on a trapping horizon, r > 2m, 

«:<1. (93) 

This means that, for a black hole of given size, the surface gravity cannot exceed a given amount. Both inequalities 
generalize those given in spherical symmetry . They may also be combined with the Penrose inequality ( pO[ ) and 
generalized to include charge |Q. 

Since {m/r)' = along a trapping horizon, ( |9l| ) implies 

to' = ■^{Airr'^y + u;(47rrV)- (94) 
an 



This is the first law of black-hole dynamics, the projection of the unified first law ( |72| ) along a trapping horizon. It 
has the same form as that of spherical symmetry , involving the radial area Airr"^ rather than the actual area A of 
a marginal surface. This suggests that dynamic black holes may in some sense admit a local Hawking temperature 



K/27r and a gravitational entropy 7rr^ |43|. It is noteworthy that, away from spherical symmetry, this entropy is 
generally not A/i. However, this is consistent with the usual result for a Kerr black hole, which has horizon area 
A = Airr"^ = 47r(f^ + o?), where r is the usual Kerr radial function [ p8[ . 

A zeroth law of black-hole dynamics, generalizing that of spherical symmetry [ ^^ , also follows from (|9l|): if r and w 
are constant on a trapping horizon, then so is the surface gravity n. Here constant r is a local equilibrium condition 
for the horizon, with constant w being a similar condition for the matter. It should be stressed that these three laws 
of black- hole dynamics all differ from the textbook laws of black- hole mechanics [|ll|,|2|, for which the zeroth and 
first laws refer to stationary black holes, specifically to a Killing horizon, while the second law refers to a physically 
unlocatable concept, event horizon. 



VIII. GRAVITATIONAL- WAVE DYNAMICS: COMPLEX POTENTIAL, WAVE EQUATION 



The shear equations (|35|), (p9[), composed into a second-order equation for the transverse conformal metric fc, 
become, taking the vacuum case, 

A(+(r2fc-i-A_)fc) =0. (95) 

This equation describes the propagation of the gravitational radiation, or the gravitational- wave dynamics. The 
original reference omitted the factor as a consequential error from (^2|). To gain more insight into this wave 
equation, it is useful to decompose the gravitational radiation into two independent polarizations. A convenient choice 
of variables turns out to be two functions (0, x) defined by 

k = e~^* sec 2x^19(81^1? + 2 tan 2x sin iJdiJ^dfp + e^* sec 2xsin^i?(i(/?(g)d(/3. (96) 

Then (f) and x respectively encode the so-called plus and cross polarizations [ p5[ of the gravitational radiation, referred 
to the quasi-spherical polar coordinates {'d,ip). Straightforward calculations show that the gravitational-radiation 
energy tensor (^ takes the form 

A„(/)A^0 + A„xA^x - i5a/35^^(A^0A5(/. + A^xA^x) 

'^aP = -. 2^ • (97) 

47r cos^ 2x 

This is an energy tensor for two coupled scalar fields (0, x)j in particular reducing to the Klein-Gordon form for 
purely plus polarization, x = 0. These two potentials may be combined into a single complex gravitational-radiation 
potential 

$ + ix- (98) 

Then the energy tensor takes the form 

47rcosh ($ - $) 



11 



where the bar denotes the complex conjugate. The overall factor involving the cross polarization $ — $ indicates 
that this is not simply the the energy tensor of a complex Klein-Gordon field <&; there is a non-trivial coupling. 
Nevertheless, the gravitational radiation has been encoded in a complex potential, acting like a complex scalar field. 
This is analogous to the gravitational potential in cylindrical symmetry, which also reduces to the Newtonian potential 
in the Newtonian limit [ p8[ . The current case differs in that $ is a potential for the gravitational radiation only, not 
including the quasi-spherical part of the geometry, for which one might instead define a gravitational potential — m/r. 
The wave equation (|9^), written explicitly in terms of {<p,x) by straightforward calculations, becomes 

A+A_0 + r"^(A+rA_0 + A+0A_r) + 2tan2x(A+xA_0 + A+<?!)A_x) = (100) 
A+A_x + r"^(A+rA_x + A+xA_r) + 2tan2x(A+xA-X - A+0A_(/i) = 0. (101) 

Written covariantly in terms of $ and the quasi-spherical wave operator 

□ = A2 + 2r"^Ar- A* (102) 

this becomes 

□ $ = 2 tanh($ - $) A$ • A»$. (103) 

This complex wave equation encodes the gravitational-wave dynamics. For purely plus polarization, $ real, the 
right-hand side vanishes and the equation is a linear wave equation, with the same form as the wave equation 
'V'^(f> = in spherical symmetry. Otherwise, the cross polarization enters the equation and consequently affects the 
wave propagation. This polarization-coupling effect is quadratic in A$ and so will be small for weak waves, thereby 
agreeing with linearized gravitational- wave theory. It should perhaps be stresse d tha t, while the gravitational-radiation 



potential $ depends on all four space-time coordinates, the wave equation ( 103 ) is an effectively two-dimensional 
equation, describing wave propagation independently at each angle of the sphere. Similarly, while the wavefronts are 
roughly spherical, the waves themselves need not be; their amplitude is independent at each angle of the sphere. 

Remarkably, the wave equation ( |103D is a type of Ernst equation, originally discovered in the quite different context 
of stationary axisymmetric space-times, with 2$ being what Ernst denoted by /i |44[ . Ernst introduced other potentials 
more commonly denoted ||l^,^,Q by 

Z = e^* (104) 

and 

£: = tanh$ (105) 



in terms of which the wave equation (102) becomes 

{Z + Z)nZ = 2AZ ■ A^Z (106) 

or 

{EE~l)nE^2EAE-A^E. (107) 

The Ernst equation has been extensively studied both in stationary axisymmetry ||l3| and in plane symmetry p6[ , 
so existing methods can be used to analyse the gravitational-wave dynamics. The current case is more complex in 
that r is generally a function of all four space-time coordinates and is not known independently of $ in the second 
approximation. However, in the first approximation, say on a Schwarzschild black-hole background, the equation 
can be treated by existing analytical methods. In particular, a direct comparison with linearized gravitational waves 
should yield insight into the non-linear effects. This at least demonstrates that, even in the very simplest case, the 
appr oximation is not some disguised version of linearized theory. In summary, the master gravitational- wave equation 
( |103| ) describes generally non-linear gravitational-wave propagation in physical circumstances where the wavefronts 
are roughly spherical. 



IX. CONFORMAL FIELDS: LOCALIZED BONDI FLUX, STRAIN 

In order to compare with observation, one needs to examine the gravitational radiation at large distances in an 
almost flat background, classically described by Bondi-Penrose asymptotic theory p^-pO|. This is straightforward in 
terms of the quasi-spherical conformal factor 
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(108) 



Then = defines conformal boundaries of the space-time, including future null infinity 3"*" and past null infinity 
for an asymptotically flat space-time. Then the conformally rescaled expansions and shears 

z9± = Q-^e± (109) 
= na± (110) 



are finite and generally non-zero at 3^. Rewriting the truncated field equations (^3[)-(|39|) of the second approximation 
yields, taking the vacuum case, 

(111) 
(112) 
(113) 
(114) 
(115) 
(116) 
(117) 

These conformal field equations are more practically suited to obtaining the gravitational waveforms at large distances, 
and are those which have been numerically implemented ||2^ ]. Similarly, the energy flux ip (fz^ ) may be rescaled to 
yield a quantity which is generally non-zero at 5^ , the conformal flux 

ifi^n-^^. (118) 

For the gravitational-radiation energy flux, the explicit expressions 



A±n = 








A±k = 




A±i9± = 












A±<;^ = 





have the same form as those for the Bondi flux at 3^ |37|. Thus the Bondi flux has been localized. The unified first 



law (^), projected along 3+, thereby reduces to a Bondi energy-loss equation. The entire preceding theory may be 
reformulated straightforwardly in terms of such conformal fields. The companion article p9| illustrates this by using 
the conformal fields, whereas this article generally uses the more familiar physical fields, such as radius r. 

The physically observable quantity to be measured by a gravitational-wave detector is the strain tensor e, which 
determines the displacements 

j^eice) (120) 

where, in Newtonian theory, e is a Cartesian basis vector in the direction of displacement. In Einstein theory, the 
displacements for a transverse vector e are given by 

h{e,e)+Sh{e,e) _ {I + Si)^ 

/i(e,e) ^2 • ^ ) 

As a transverse surface approaches a metric sphere at 5^, Sfl/fl and fc(e, e) 1, so the formula yields 

^<5fc(e,e)^2. (122) 

This leads to the identification of e with 6k/2, or more precisely with 

i f LLakdT^]- [ Vl{a+ ^+ + a- ',-)dT (123) 
^ J ^ ^ J J 

where 7 is a worldline normal to the transverse surfaces and a = d/dr = a+Z+ -I- a~L is a vector tangent to 7. Since 
the strain itself vanishes at 5^, it is more practical to introduce a quantity which is finite at future null infinity 
the conformal strain tensor 
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e^l U- dx- (124) 
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2 

where 7 is now a null curve of constant (x"*", "i?, (/?). Then 

e = - (125) 
r 

is the strain tensor at a large distance r from the source. Thus the variables of the quasi-spherical approximation 
are directly related to the observable strain. There is no need to match with another far-zone approximation, such 
as a linearized approximation. For instance, graphs of the strain waveforms, e against x~ at a given angle 93), are 
produced as output by the existing numerical codes [ p8| . 

A familiar asymptotic relation between strain and energy density is recovered: the energy density of outgoing 
gravitational radiation (<j_|_ = 0) near S"*" is, taking the normalization = l/v^, 

= (126) 

where one should remember that the two independent components of e each contribute twice in the norm. Regarding 
conventions, an unfortunate relic of linearized theory is a definition of amplitude of a gravitational wave as that 
of the relevant component of the transverse traceless metric perturbation [psf, which would correspond here to k. 
As above, this differs by a factor of two from the strain amplitude, which is a much more natural definition of 
gravitational-radiation amplitude, being the observable quantity. This is also simply related to the amplitude of the 
gravitational-radiation potential $ by 

A $A $ 

e(a, a) = ' , ' (127) 

^ ' 87rcosh2($- $) ^ ' 

and is equal for purely plus polarization, $ = $, as for Einstein- Rosen gravitational waves |Q. Another reason for 
the latter convention is that the energy density has the same form and numerical factor, (amplitude x frequency) ^/Stt, 
as an electromagnetic wave in natural units. 



X. CONCLUSION 



This article has described an astrophysically realistic approximation scheme in which both gravitational radiation 
and black holes are locally defined, along with their physical attributes, with each dynamically influencing the other. 
The approximation holds where the gravitational wavefronts and black hole (or other astrophysical object) are roughly 
spherical. For black holes, this is an extension of the general theory of black-hole dynamics |^-|3^ and its previous 
formulations in spherical | |3^ and cylindrical symmetry psf . In particular, local definitions of gravitational mass- 
energy m ( ^8| ) and surface gravity k (|70| ) have been given, satisfying zeroth, first ( p^ ) and second (|87| ) laws of 
black-hole dynamics. For gravitational radiation, this is a new kind of approximation, quite different from the 
commonly used asymptotic and high-frequency (or merely linearized) approximations. Most crucially, gravitational 
radiation is localized in this approximation, something which is usually argued to be impossible in general and in other 
existing approximations, unless one counts symmetric space-times [p8| . The general non-localizability of gravitational 
waves has long been an obstacle to understanding the physics of their most interesting, dynamic, strong-gravity 
sources. Remarkably, in the quasi-spherical approximation, the gravitational-wave dynamics is described simply by 
an effectively two-dimensional wave equation (103). This analytically tractable Ernst equation can be used to study 
non-linear wave-propagation effects, for instance close to a black hole. 

A local energy tensor O (43) for the gravitational radiation has been defined, taking a scalar-field form in terms of 
a complex gravitational-radiation potential (|9^). Including Q like a matter energy tensor in the truncated Einstein 
equations thereby describes gravitational radiation reaction, the back-reaction of the radiation on the space-time, e.g. 
on the black ho le p roducing it. A reliability test is provided by comparing this second approximation with the first 
approximation ||27| , in which radiation reaction is not included and the radiation propagates on a background. In the 
second approximation, the gravitational-radiation energy tensor Q enters a covariant conservation law, conservation 
of energy (|8^). Specifically, with respect to a canonical flow of time, the combined energy- momentum of the matter 
and gravitational radiation is a Noether current, with the Noether charge being the mass-energy m. The energy 
conservation law can be expressed as a local first law ( [72[ ) equating changes in m to work and energy-supply terms, 
including the energy flux ( |7^ ) of the gravitational radiation. This unified first law includes both the first law of 
black-hole dynamics and a Bondi energy-loss equation. 
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It may be useful to summarize here how the Einstein equation has been simphfied and physically understood. Of the 
ten components of the full Einstein equation, four components of the truncated equations appear to involve angular 
mom entum but are irrelevant to the radiation problem. Two of the remaining components constitute the wave equation 



( 103 ) describing the propagation of the gravitational radiation. The remaining four components are the quasi-spherical 
equations, of which one is an integrability condition for the others, at least in vacuo and for some matter fields, such 
as a Klein-Gordon field. Of these, one is a relativistic version of Newton's inverse-square law of gravitation, (|9l|), 
and two are equivalent, up to an integration constant, to the first law ([72|). Thus a comprehensive physical picture 
has emerged: the Einstein equation has been essentially reduced to a quasi-Newtonian law of gravitation, energy 
conservation and a propagation equation for gravitational radiation. 

As an example of current interest to many researchers, the approximation can be used to describe gravitational 
radiation from a roughly spherical but dynamically evolving black hole formed by binary black-hole coalescence. The 
distorted black hole emits gravitational radiation, absorbs backscattered radiation, thereby increases in area and 
traps some outgoing radiation, changes shape accordingly and consequently emits more radiation, and so on. This 
ongoing dynamical process may now be described in a local way by physically interpretable equations, involving 
physical quantities such as black-hole mass and surface gravity and gravitational-radiation energy flux. Apart from 
such practical applications to gravitational-wave astronomy, this provides a rich arena in which to advance physical 
understanding of both gravitational radiation and black holes, and the dynamical interaction between them. 
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